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PACS 74 . 70 . -b - Superconducting materials 

PACS 74 . 25 . Jb - Electronic structure 

PACS 74 . 25 . Ha - Magneticproperties 

PACS 71 .20.-b - Electron density of states and band structure of crystalline solids 

Abstract. - By first-principles calculations, we present a doping-dependent phase diagram of 
LaOMAs (M=V-Cu) family. It is characterized as antiferromagnetic semiconductor around 
LaOMnAs side and ferromagnetic metal around LaOCoAs. Both LaOFeAs and LaONiAs, where 
superconductivity were discovered, are located at the borderline of magnetic phases. Extensive 
Fermi surface analysis suggests that the observed superconductivity is of electron-type in its ori- 
gin. We discuss possible pairing mechanisms in the context of competing ferromagnetic phases 
found in this work and the ferromagnetic spin fluctuations. 
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The studies on new superconductors, particularly non- 
cupcrate layered compounds, are always exciting and open 
up new windows for the possible further raising of transi- 
tion temperature T c . Except few of non-transition- metal 
compounds, such as MgB2, where superconductivity with 
T c up to 39K was found [1], problems in many cases are 
as complicated and challenging as what we found in cu- 
pcratcs. In the layered ruthenates, the ferromagnetic spin 
fluctuation is important and spin-triplet p-wave charac- 
ter was suggested [2] for the superconductivity found in 
Sr2Ru04 with T c ~1K [3]. In a more recent example 
Na 2; Co02-(H20)j / [4], the geometry fluctuations due to 
triangle-lattice are extensively discussed [5]. Here we will 
show that a rich doping-dependent phase diagram can be 
realized in the new family of layered compounds LaOMP 
or LaOMAs (M=V-Cu), where up to T C =26K supercon- 
ductivity was reached very recently in LaOFeAs after F~ 
doping [6-8]. 

The quaternary oxypnictides LaOMAs crystallize in lay- 
ered tetragonal structure with PA/nmm symmetry [9]. 
Each transition-metal (oxygen) layer is sandwiched by two 
nearest-neighbor As (La) atomic layers, which form edge- 
shared tetrahedrons around the M (oxygen) sites. The 
(MAs)~ and (LaO) + triple-layer-subgroups stack alterna- 
tively along the c-axis. The positions of La or As sheets 
are determined by two internal parameters, zj, a and za s . 



which define the inter-layer distances of La-0 and M- 
As, respectively. It is important that this series of com- 
pounds are chemically stable such that systematical tun- 
ing is available without altering the structure and sym- 
metry significantly. For instance, a variety of compounds 
can be synthesized by the replacement of transition-metal 
elements, where both the electron doping and hole doping 
can be realized by replacing 2_ or La 3+ ions. Except 
the early report for the structure study [9], the detailed 
studies on the electronic and magnetic properties for this 
series of compounds are still in its early stage. It was first 
reported in 2006 that superconductivity can be realized in 
LaOFeP below 4K, and the T c was increased to 7K by F- 
doping [6]. Later, superconductivity with T c about 2K was 
reported for LaONiP [7], and T c around 26K was reached 
very recently in LaOFeAs again after F-doping [8] . We will 
present in this letter that both M=Fe and Ni compounds 
locate at special positions of the global phase diagram for 
the series of M substituted compounds. The competing 
magnetic and superconducting phases found in our global 
phase diagram provide important clues on the possible 
pairing mechanism in this class of materials. 

The phase diagram is constructed from first-principles 
calculations based on density functional theory with gen- 
eralized gradient approximation (GGA) of PBE-type [10] 
for the exchange-correlation potential. We use the plane- 
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Table 1: The optimized lattice parameters (a and c), internal 
coordinates (zLa and zas), and calculated specific heat coef- 
ficient 70 and bare susceptibility xo f° r series of compounds 
LaOMAs in the NM state. 





o(A) 


c(A) 


ZLa 


ZAs 




Xo(10- 5 S) 


V 


3.9965 


9.3655 


0.1367 


0.1643 


7.47 


10.25 


Cr 


3.9355 


9.3467 


0.1397 


0.1592 


9.06 


12.43 


Mn 


4.0355 


8.7822 


0.1438 


0.1427 


11.37 


15.59 


Fe 


4.0396 


8.6289 


0.1461 


0.1369 


5.52 


7.57 


Co 


4.0621 


8.5434 


0.1466 


0.1369 


7.41 


10.17 


Ni 


4.1404 


8.3139 


0.1466 


0.1366 


3.81 


5.23 


Cu 


4.1442 


8.5681 


0.1424 


0.1539 


3.98 


5.46 



wave pseudopotcntial method, and the ultra-soft pseu- 
dopcntentail scheme [11] is adopted. The convergence of 
total energy calculations with respect to number of K- 
points and cut-off energy (of plane wave expansion) is well 
checked, and final results are double checked using the 
full-potential linearized augmented plane-wave (FLAPW) 
method (WIN2K package) [12]. The series of LaOMAs 
compounds with M ranging from V to Cu are all stud- 
ied with full lattice optimizations using experimental 2 
Fe cell, and the non-magnetic (NM) , ferromagnetic (FM) , 
and (intra-layer) antiferromagnetic states are treated. The 
same approach has been also applied to LaOMP series, 
qualitatively same results are obtained, we therefore con- 
centrate our following discussions on LaOMAs series for 
consistence. 

As shown in Table I, the optimized lattice parame- 
ters and internal coordinates for Fc and Ni compounds 
are in excellent agreement with available experimental 
data [6-9] , which demonstrates the quality of our present 
calculations. The optimized parameters are used in our 
calculations for all the compounds. Moving from V to Cu, 
the lattice parameters arc only slightly modified, despite 
of the dramatical change of number of 3d electrons (from 
d 3 of V 2+ to d 9 of Cu 2+ ), suggesting the less sensitivity of 
lattice distortion. 

Fig.l shows the phase diagram computed for the whole 
range of compounds. The solid lines and the dashed lines 
represent the stabilization energies of FM and AF states 
relative to NM solution, respectively. Two distinct phase 
regions can be identified: the dome below Fe gives AF 
ground state, while the one around Co is ferromagnetic. 
The computed magnetic moments (see Fig.l) show that 
the left-hand side AF phase region are strongly spin po- 
larized, while the magnetic moment of the ferromagnetic 
phase around Co is small, suggesting possibly different ori- 
gins of two magnetic phase regions. It is important to note 
that both LaOFeAs and LaONiAs arc located at special 
positions of the phase diagram: Fe compound is at the 
borderline between the AF and FM phase regions; while 
Ni compound is at the other border of the FM phase. 
The overall phase diagram clearly suggests that the mag- 
netic instabilities are crucially important to understand 
the physical properties of this new family of compounds. 
The present calcualtions are done based on experimental 
crystal cell with 2 Fe per cell [8]. However, we notice that 
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Fig. 1: Phase diagram: (a) The calculated stabilization ener- 
gies of FM and AF states with respect to NM state for various 
compounds with increasing number of 3d electrons, (b) The 
corresponding magnetic moments obtained from calculations. 

strong nesting effect exists particularly for LaOFeAs. This 
will lead to stripe-type spin-density-wave (SDW) ground 
state with \/2 x \/2 super-cell structure. The detailed re- 
sults for this SDW state will be presented in a separate 
paper combined with experimental results [13]. 

To understand the electronic structures, we show the to- 
tal and projected density of states (DOS) of various com- 
pounds in Fig. 2. Let's start from LaOFeAs (Fig. 2(a)). 
The states between -2eV and +2eV are mostly from Fe- 
3d states, just below which are the states of O-p and As-p 
(from -6eV to -2eV). The p-d hybridization between O and 
Fe is negligible, while that between As and Fe is sizable. 
This As-Mp-d hybridization is enhanced by changing to Ni 
compound (see the projected DOS of Fig. 2(b)). Since the 
transition metal sites are coordinates by As-tetrahedron, 
the crystal field will normally splits the five d orbitals into 
low-lying two-folds e g states and up-lying three-folds t 2g 
states. However, the As-tetrahedron are actually much 
distorted from its normal shape (squeezed along c by about 
20%). This distortion will further split the e g and t 2g man- 
ifolds significantly making the final orbital distributions 
complicated. As the results, what actually happens is op- 
posite to what we expect from simple tetrahedra crystal 
field: the low-lying manifold is three folds and the higher- 
lying manifold is two folds, between which a pseudo-gap 
about 0.5eV exists. For Fe compound, the nominal num- 
ber of d electrons is 6, and the low- lying three- fold bands 
are nearly fully occupied with the Fermi level Ef located 
very close to the deep of the pseudo-gap. If the crys- 
tal field is strong enough, which may be achieved by the 
substitution of As atoms, a simple band-insulator will be 
expected by enlarging the pseudo-gap. From the calcu- 
lated DOS at Ef, we estimate the bare susceptibility and 
specific heat coefficient, which are xo=7-6xlO -5 emu/mol 
and 7o=5.5 mJ/K 2 mol for Fe compound (see Table I for 
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other compound). Taking the experimental susceptibility 
about 50xl0 -5 emu/mol at 300K, a renormalization fac- 
tor about 6.6 is suggested, which is close to that shown in 
Sr 2 Ru0 4 [14]. 

Moving to Mn compound, the half filled d-shell (about 
d ) will gain energy from the Hund's coupling, and the 
spin polarized states will be favored as the results. In re- 
ality, the calculated total energies of both FM and AF 
solutions are much lower than that of NM solution (about 
0.3eV/f.u for the former and 0.55eV for the later). The 
reason why the AF solution is more favored is that a gap 
about 0.2eV is opened in the AF solution (see Fig. 2(e)). 
It is mentioned (but without showing data) in the recent 
experimental paper [6] that LaOMnP is a semiconduc- 
tor. This is consistent with our prediction, and further 
we show that the ground state of Mn compound is AF. 
The calculated spin moment is about 3.1/is/Mn, which 
is much reduced from the expected 5[Ab /Mn of high spin 
state. There are two possibilities to explain the reduced 
moment. One is that Mn is in intermediate spin state 
rather than high-spin state. If this is the case, certain 
kinds of orbital ordering would be expected for an AF in- 
sulator. However, our calculated occupation numbers of 
projected 3d-shell orbitals are quite uniform, suggesting 
this possibility is unlikely. The second possibility is due 
to either p-d or d-d hybridization (particularly for nar- 
row gap system). The calculated spin moment is about 
4.3/is/Mn for typical high-spin AF insulator MnO [15], 
where only the p-d hybridization is important. However, 
here we point out that the Mn-Mn distance in LaOMnAs 
is about 2.8A, which is much shorter than what was found 
in MnO (about 3.2A), and is actually very close to the dis- 
tance in elementary-Mn (about 2.6~2.7A). The direct d-d 
hybridization will be much enhanced by such short dis- 
tance, which will again reduce the moment. The effect of 
direct d-d overlap has been addressed in previous study for 
LaOFeP [16], however we emphasize here that this is gen- 
erally true for all the compounds of this family as shown 
by the optimized structure (Table I). 

The stabilization of FM phase region at the right hand 
side has different origin as will be discussed here for LaO- 
CoAs. Co has one more d electron than Fe, therefore the 
Fermi level is lifted up and located above the pseudo-gap. 
What is interesting is that a strong Van-Hove singular- 
ity (VHS) is present just at the Fermi level of LaOCoAS 
NM DOS as shown in Fig.2(d). The high N(E f ) in the 
presence of this VHS will push the system to be itinerant 
FM due to Stoner instability. This mechanism is further 
supported by the following factors: (1) the FM region is 
relatively narrow; (2) the polarized spin moment is small 
(about 0.5/zb/Co); (3) the energy gain is also small (about 
35mcV/f.u.). By adding one more electron, for LaONiAs, 
the Fermi level is shifted away from the VHS, the system is 
recovered to be NM again as shown in the phase diagram. 

Having finished the discussions for the phase diagram 
and the general picture of the electronic structure, now let 
us focus on the LaOFeAs and LaONiAs, where supercon- 
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Fig. 2: The calculated electron density of states (DOS) for 
various compounds in different states (see text part for expla- 
nation) . 



ductivity are discovered. Fig. 3 gives the calculated band 
structures and Fermi surfaces (FS) for both compounds 
(the window for d-bands is shown, and there are totally ten 
d bands in our unit cell). First of all, the band dispersions 
along the z direction are all very weak suggesting the 2- 
dimcnsional nature of those compounds. Considering the 
in-plane dispersions, the band structure of LaOFeAs can 
be schematically separated into two parts. The bands be- 
low Ef are relative flat and have little contribution com- 
ing from the As-p states, while the bands above Ef are 
quite dispersive (except some flat branches around +leV 
which corresponds to the VHS discussed above), and have 
large weight coming from the As-p character (as shown 
by the projected fat-bands plot). Those dispersive bands 
form electron-like FS cylinders around the M-A lines of 
the Brillouin zone, and hole-like FS cylinders are formed 
around the T-Z lines due to the Fermi level crossing of the 
bands from the lower part. For LaONiAs, all the flat bands 
are pushed down to below Fermi level, and only the dis- 
persive bands remain to cross the Fermi level, which give 
the large electron FS around M-A lines, but hole-type FS 
around X-R line instead of T-Z. 

The experimental results show that the LaOFeAs is su- 
perconducting only after electron- type F-doping [8], and 
also the T c of LaOFeP (which has very similar band struc- 
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Fig. 3: The calculated band structures of (a) LaOFeAs and (b) 
LaONiAs. The fat-bands are shown with projection to As-p 
states. The lower panels show the calculated Fermi surfaces of 
LaOFeAs after 10% F-doping (left) and LaONiAs (right). 



ture [16] as LaOFeAs) is enhanced with F-doping [6]. All 
those evidences suggest that the electron-like FS formed 
by the dispersive band should be responsible for the su- 
perconductivity. Taking the calculated total DOS of 
LaOFeAs, the N(Ef) will decrease with up-shifting of 
Fermi level, in opposite to the experimental trend observed 
for F-doping. However, if we only take into account the 
contribution from the electron-like FS, its DOS will in- 
crease with increasing Fermi energy. Actually, as shown 
in Fig.3, hole-like FS of LaOFeAs after 10% F-doping are 
much reduced, and electron-like FS are enlarged, com- 
pared to the case without doping [16]. The very recent 
Hall measurements also suggested the electron-type con- 
ductivity in LaONiAs after F-doping [17]. 

Up to now, superconductivity has been discovered in 
several transition metal systems, including the cuperates, 
ruthenates and cobaltates. Most of these superconducting 
materials are hole doped and only a few of them are elec- 
tron doped. In cuperates, the electron doped compounds 
usually have much lower T c than the hole doped ones. In 
sodium cobaltates, only hole-doped compounds are real- 
ized. The present systems are also layered compounds, 
and the p-d hybridization is sizable. Considering the close 
vicinity to the ferromagnetic instability, the present sys- 
tems are very similar to Sr2Ru04, where strong ferromag- 
netic fluctuations favor triplet pairing. In such case, the 
question is whether the pairing state is unitary and time 
reversal invariant, analogous to the B phase of 3 He, or 
non-unitary and time reversal breaking, analogous to the 



A phase of 3 He [18]. Our calculations shows that LaOFeAs 
can be basically characterized as low density electron car- 
riers doped on top of a band insulator with filled d 6 valence 
orbitals. In 3 He, the B phase is realized under the low or 
ambient pressure condition while A phase is realized under 
high pressure condition close to solidification. By analogy, 
we suggest that the low electron density system LaOFeAs 
is in the weak coupling limit, which generally favors the 
unitary, or B phase like, pairing symmetry. Upon fur- 
ther increasing doping, the non-unitary A phase could be 
realized. In two dimensions, unitary B state can be char- 
acterized as a state where the up (down) spin electrons 
are paired in the p x + ip y (p x — ip y ) state, so that the 
time reversal symmetry is preserved. This state is similar 
to the topologically non-trivial state, characterized by Z2 
invariant [19], found in quantized spin-Hall insulator [20]. 
In contrast, the time reversal symmetry breaking A phase 
could be realized in Sr2Ru04 partly because the carriers 
density is high. 

The topological nature of the proposed pairing state 
for LaOFeAs implies the existence of counter-propagating 
edge states which can be tested experimentally. In the 
bulk, the pairing state is fully gapped, and STM exper- 
iment would show a full gap in the I-V characteristics. 
However, moving to the edge, the STM experiment would 
show a gapless spectrum, revealing the gapless edge states 
protected by the time reversal symmetry. Most strikingly, 
in the presence of a magnetic impurity near the edge, the 
local density of states would show a gap again, due to 
breaking of time reversal symmetry 
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